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SUMMARY 


The  shock-on-shock  phenomenon  represents  a potentially  dangerous 
environment  to  high  performance  strategic  interceptors.  Therefore,  there 
is  a need  to  have  the  capability  to  predict  the  pressure  distribution 
history  due  to  such  encounters.  The  objective  of  this  program  was  to 
develop  a computerized  model  which  is  inexpensive  to  operate  as  well  as 
accurate.  The  knowledge  required  to  simplify  the  model  is  derived  from 
the  Defense  Nuclear  Agency  sled  test  data  and  the  NASA  Ames  rigorous 
analytical  results. 

The  model  presented  here  permits  prediction  of  the  circumferential 
distribution  of  peak  pressures  that  result  from  an  encounter  at  any 
angle  from  nose-on  to  broadside.  The  sharp  supersonic  conical  body  may 
also  be  at  a small  (less  than  the  cone  half-angle)  angle  of  attack.  The 
model  also  allows  prediction  of  the  circumferential  distribution  of  the 
duration  of  the  shock-on-shock  pulse  which  is  defined  as  the  time  it 
takes  for  the  surface  pressure  to  reach  the  "quasi-steady"  pressure 
level  associated  with  complete  engulfment  of  the  vehicle.  Finally,  the 
model  permits  estimating  the  engulfment  time  of  a particular  axial  station 
of  the  vehicle,  i.e.,  the  time  interval  between  the  arrival  of  the  wave 
front  at  the  subject  station  on  the  windward  side  of  the  cone  and  the 
arrival  at  that  station  on  the  leeward  side. 

Comparisons  between  the  new  model  and  the  existing  data  show  good 
agreement  and,  where  no  comparisons  are  possible,  the  new  results  are 
plausible.  In  spite  of  the  simplifying  assumptions  and  some  empiricism, 
the  model  is  expected  to  be  effective  over  a fairly  large  range  of  blast 
and  vehicle  conditions  except  for  angle  of  attack.  When  the  angle  of 
attack  increases  to  the  point  of  flow  separation  on  the  leeside  of  the 
cone,  the  model  no  longer  applies;  it  could,  however,  be  modified  to  do  so. 
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1.0  INTRODUCTION 


The  shock-shock  interaction  phenomenon  occurs  when  a supersonic 
vehicle  encounters  a shock  wave.  Recent  interest  in  this  problem  is 
due  to  the  high  probability  of  such  an  occurrence  during  the  defense 
of  the  Minuteman  system  by  an  interceptor.  The  threat  of  many  reentry 
vehicles  attacking  a single  target  necessarily  results  in  a fratricidal 
environment  being  created  for  an  interceptor  by  the  detonations  of  pre- 
ceding interceptors.  Such  an  encounter  results  in  high  frequency  shocks 
transmitted  through  substructure  which  can  damage  components.  The  sur- 
face pressure  is  usually  maximum  immediately  behind  the  wave  which 
strikes  the  surface  of  the  vehicle  and  then  it  decays  to  a "quasi- 
steady" value.  The  latter  pressure  is  referred  to  as  quasi-steady 
because  the  apparently  steady  conditions  will  actually  vary  as  the 
vehicle  traverses  the  blast  sphere.  However,  relative  to  the  short 
time  it  takes  for  the  vehicle  to  be  engulfed  by  the  wave,  the  conditions 
behind  the  blast  front  may  be  assumed  to  be  steady  and  the  surface 
pressure  similarly  achieves  a steady  value.  The  time  to  decay  to 
quasi-steady  pressure  is  defined  as  the  pulse  duration. 

The  objective  of  this  effort  was  to  produce  a relatively  simple, 

but  accurate  model  of  the  shock-on-shock  phenomenon  as  it  occurs  when 
a sharp  supersonic  conical  vehicle  encounters  a planar  blast  wave. 

The  desire  is  to  be  able  to  predict  the  following  features  of  the 
phenomenon,  namely,  the  peak  pressures  and  pulse  durations  produced 
along  the  windward  and  leeward  rays  of  the  cone  and  the  circumferen- 
tial distributions  of  these  quantities.  Any  encounter  angle  from 
nose-on  to  broadside  is  permitted  and  the  cone  may  be  at  small  (less 
than  the  cone  half-angle)  angles  of  attack.  The  requirement  for 
simplicity  stems  from  the  need  for  a low  cost  computer  model. 
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The  approach  taken  in  structuring  the  model  was  to  formulate  a set 
of  simplified  assumptions  based  on  the  knowledge  gained  from:  1)  the 

rigorous  analytical  results  of  Reference  1,  2)  the  test  results  of 
Reference  2,  and  3)  Reference  3 which  presents  a comparative  study 
between  these  sources  and  established  the  framework  for  the  model  pre- 
sented herein. 

The  following  description  of  the  problem  (Section  2.0)  is  excerpted 
from  Reference  3.  Section  3.0  summarizes  the  current  model  and  pre- 
sents comparisons  with  other  data  sources.  A listing  is  presented  in 
Section  4.0  along  with  user  information  pertaining  to  the  computerized 
version  of  the  program. 
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2.0  THE  PROBLEM  - A FUNCTION  OF  ENCOUNTER  ANGLE 

The  encounter  between  a blast  wave  and  a supersonic  vehicle 
produces  a highly  transient  pressure  pulse  which  varies  in  strength 
and  duration  depending  on  a combination  of  parameters:  geometry  and 

flight  condition  of  the  vehicle  and  the  strength  and  orientation  of  the 
blast  wave  relative  to  the  surface  of  the  vehicle.  Interest  here  is 
limited  to  sharp  cones  so  the  vehicle  condition  at  a given  altitude 
is  defined  by  a combination  of  cone  angle  (6C) , angle  of  attack  (a) , 
and  vehicle  velocity  (V) . The  strength  of  the  blast  wave  is  defined 
by  the  ratio  of  the  pressure  behind  to  the  pressure  in  front  of  the 
wave  (P2/P1)  as  it  travels  through  the  atmosphere.  The  encounter 
angle  (6)  is  defined  as  the  angle  between  the  blast  wave  and  the 
perpendicular  to  the  vehicle  axis  as  shown  in  Figure  1. 


Region  (?) 


Figure  1.  Blast  Wave  - Vehicle  Encounter  Geometry 
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The  primary  determinant  in  the  level  of  pressure  is  the  strength 
of  the  blast  wave,  but  the  encounter  angle  (3)  also  plays  an  important 
role.  As  6 varies  from  nose-on  to  broadside,  several  distinct  classes 
of  encounters  can  be  identified.  One  subdivision  of  encounters  is 
marked  by  the  bow  wave  angle.  For  8's  between  nose-on  (8=0)  and 
parallel  to  the  bow  wave  angle  (i.e.,  8=90  - 0W) , the  intersection 
of  the  blast  and  bow  waves  during  engulf ment  moves  from  the  apex 
of  the  cone  aftward.  For  8's  between  (90  - 6W)  and  (90  - 0C) , the 
blast-bow  wave  intersection  moves  forward  along  the  bow  wave.  The  two 
cases  must  be,  in  general,  treated  separately. 

The  range  of  8's  between  nose-on  and  (90  - 0W)  must  be  subdivided 
further  because  the  nature  of  the  wave  pattern  which  develops  on  the 
surface  varies  from  Mach  reflection  to  regular  reflection  as  the 
encounter  varies  from  nose-on  towards  broadside.  The  two  types  of 
wave  pattern  are  illustrated  in  Figures  2a  and  2b  respectively. 

Unless  a solution  of  the  complete  flow  field  is  developed,  as  was 
done  in  Reference  1,  the  two  types  of  wave  patterns  must  be  treated 
separately  because  there  is,  otherwise,  no  precise  way  of  predicting 
the  location  and/or  height  of  the  Mach  stem.  The  regular  reflection 
pattern,  on  the  other  hand,  can  be  solved  rigorously  as  long  as  the 
transmitted  and  reflected  waves  are  straight,  i.e.,  immersed  in 
uniform  flow  fields.  As  discussed  in  Section  3.1,  the  presence  of 
gradients  in  the  flow  field  leads  to  bending  of  the  waves  which 
complicates  the  solution  of  the  wave  pattern.  It  is  shown  that 
acceptable  approximations  are  possible,  however,  and  that  good  agree- 
ment is  achieved  between  the  simple  model  and  the  rigorous  solution 
of  Reference  1. 

Encounters  near  broadside,  where  (90  - 0W)  <_  8 _ (90  - 0C) , 
must  also  be  subdivided  into  two  subclasses  depending  on  the  orien- 
tation of  the  waves  which  emerge  from  the  blast  bow  wave  interaction. 
While  the  blast-bow  wave  intersection  moves  forward,  the  point  of 


8 


intercept  at  the  vehicle  surface  may  move  from  front  to  back  depending 
on  the  orientation  of  the  wave  transmitted  to  the  surface  as  shown 
in  Figure  3. 


A regular  reflection  will  occur  in  either  case,  but  any  attempt  to 
model  the  phenomena  must  reflect  the  differences  to  correctly  describe 
the  velocity  of  the  transmitted  wave  relative  to  the  gas  at  the  vehicle 
surface. 

The  complexities  associated  with  the  various  subclasses  of  en- 
counters are  manifested  in  the  development  of  models  of  the  shock- 
shock  phenomena.  The  earliest  analytical  success  occurred  for  nose-on 
encounters  which  are  further  divided  into  two  subclasses:  sharp  and 
blunt  vehicles.  The  nose-on  encounter  of  a slender  sharp  cone  was 
treated  by  Inger  in  Reference  4 and  the  blunt  body  was  treated 
successfully  by  McNamara  and  Taylor  in  References  5 and  6,  respectively. 
It  was  possible  to  demonstrate  the  success  of  these  solutions  because 
of  the  test  data  available  (e.g..  References  7 and  8).  Meanwhile, 
solutions  were  also  developed  for  oblique  encounters,  but  remained 
unverifiable  until  the  recent  sled  test  data  became  available. 

Probably  the  highest  quality  experiments  involving  oblique  encoun- 
ters produced  previously  were  the  White  Oak  Laboratory  of  the  Naval 
Surface  Weapons  Center  using  a ballistic  range  in  combination  with  a 
shock  tube.  Excellent  photographic  data  (References  9 and  10)  were 
taken,  but  the  use  of  ballistic  models  precluded  the  acquisition  of 
actual  pressure  measurements.  Thus,  actual  pressure  measurements  were 
limited  essentially  to  nose-on  encounters  because  of  the  limitations 
of  the  test  techniques. 

Interestingly,  the  sled  tests  cannot  be  used  to  simulate  nose-on 
encounters  because  of  interference  between  the  wave  and  the  ground. 
However,  the  sled  test  scheme  does  permit  simulation  of  Mach  stem, 
regular,  and  broadside  encounters,  thereby  filling  a large  gap  in  the 
body  of  test  data. 


3.0  SHOCK-ON-SHOCK  MODEL 


The  description  of  the  model  discussed  in  this  section  is  sub- 
divided in  segments.  The  peak  pressure  on  the  windward  ray  of  the  cone 
is  analyzed,  followed  by  the  model  of  the  peak  pressure  on  the  leeward 
ray.  The  corresponding  analyses  of  the  windward  and  leeward  pulse 
durations  are  then  presented,  followed  by  a discussion  of  the  circum- 
ferential distributions. 

The  rationale  presented  here  is  used  for  both  zero  and  non-zero 
angle  of  attack  (o)  by  replacing  the  real  cone  angle  (0C)  by  an  effec- 
tive cone  (Figure  4).  The  windward  ray  analysis  was  performed  on  an 
effective  cone  with  semi-vertex  angle  of  0c+a  at  zero  angle  of  attack, 
while  in  the  leeward  ray  analysis,  the  semi-vertex  angle  is  0c-a. 
Similarly  the  bow  wave  angle  is  measured  from  the  centerline  of  the 
effective  cone. 
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3.1  Peak  Windward  Pressure 


The  starting  point  for  the  analysis  was  to  examine  some  results 
taken  from  Kutler  and  Sakell  (Reference  1)  and  shown  in  Figure  5. 
Figure  5 shows  the  isobar  maps  calculated  with  the  Kutler  solution 
for  an  11.2-degree  half-angle  cone  at  Mach  5.  The  isobars  show  the 
shock  wave  patterns  that  result  when  the  cone  encounters  a blast  wave 
of  strength  1.6  (i.e.,  at  six  different  encounter  angles.  At 

the  nose-on  encounter  (B=0),  the  wave  transmitted  to  the  surface  is 
very  nearly  straight.  As  the  encounter  angle  becomes  oblique,  the 
transmitted  wave  is  seen  to  be  curved.  It  was  decided  in  formulating 
the  new  model  referred  to  here  that  the  curvature  in  some  form  should 
be  considered. 


Figure  5.  Windward  and  Leeward  Isobar  Maps  at  Selected  Encounter  Angles 
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Careful  examination  of  the  isobar  maps  shows  that  the  height  of 
the  stem  formed  on  the  windward  ray  of  the  cone  is  relatively  small 
in  comparison  to  the  shock  layer  thickness.  The  reflected  wave 
bends  away  very  quickly  from  being  perpendicular  to  the  surface  to 
assume  an  orientation  which  is  compatible  with  the  new  bow  wave  forming 
on  the  vehicle  as  a result  of  being  immersed  in  the  post-blast  flow 
field.  In  fact,  at  all  of  the  calculated  cases  except  nose-on,  the 
reflected  wave  bends  almost  immediately  upon  leaving  the  surface. 

At  some  encounter  angle  between  nose-on  and  40  degrees,  the 
wave  pattern  changes  from  Mach  stem  to  regular  reflection.  Further- 
more, as  might  be  expected,  there  is  no  abrupt  change  in  the  wave 
pattern  as  the  encounter  angle  changes  from  one  which  produces  regular 
reflection  to  one  which  produces  a Mach  stem.  Stated  another  way,  it 
appears  that  the  transmitted  wave  shape  and  orientation  right  down 
to  the  surface  may  be  determined  without  regard  for  the  type  of 
reflection  which  will  occur.  If  the  transmitted  wave  is  described 
properly,  then  the  test  for  regular  reflection  is  made  at  the  surface 
and  the  peak  pressure  may  be  calculated  accordingly,  i.e.,  as  the 
pressure  behind  a regular  reflection  or  as  the  pressure  behind  a small 
normal  shock  located  at  the  point  of  impingement  of  the  transmitted 
wave  on  the  surface. 

The  queston  is  how  to  solve  for  the  bent  transmitted  wave  in  a 
simple  fashion.  The  Kutler  solution  accomplishes  this  by  subdividing 
the  shock  layer  into  many  small  segments  and  solves  for  the  changes 
in  flow  conditions  from  point  to  point  in  the  grid.  The  many  points 
in  the  grid  result  in  a more  rigorous  solution  and  increased  computer 
cost.  Noting  that  the  transmitted  wave  actually  is  not  highly  curved, 
it  was  decided  that  the  curved  wave  might  be  reasonably  approximated 
with  two  straight  segments  as  indicated  in  Figure  6.  An  iterative 
procedure  is  used  to  solve  for  the  outer  segment  of  the  transmitted 
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wave  by  solving  the  two-shock  interaction  at  the  point  where  the  bow 
and  blast  waves  intersect. 

Both  waves  are  assumed  to  be  locally  two-dimensional  and  their 
interaction  results  in  two  waves  which  are  the  beginnings  of  the  post- 
encounter bow  wave  and  the  transmitted  wave,  respectively.  The 
rationale  used  in  solving  the  interaction  of  two  planar  waves  is: 

]L  Assume  a value  for  pressure  jump  across  wave  P-R  (Figure  6) 

1_  Calculate  the  flow  direction  in  region  4 (Figure  6) 

_3  Calculate  the  flow  direction  in  region  5,  assuming  the 
pressure  in  region  4 

If  the  flow  directions  in  regions  4 and  5 (Figure  6) 
agree  then  the  solution  is  found,  otherwise  the  procedure 
is  repeated  with  a new  pressure  jump  assumed  in  step  one. 
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The  key  equations  defining  the  above  procedure  are  summarized 
below.  To  solve  for  the  system  of  shock  waves  which  result  from  the 
intersection  of  the  blast  and  bow  waves,  it  is  useful  to  place  the 
frame  of  reference  at  the  point  of  intersection  (P,  Figure  6).  The 
velocity  of  point  P,  along  the  bow  wave,  relative  to  the  body  fixed 
coordinates  (^,yw>  is 


VP/W 


cos  (g-6) 
vB/Wt  cos  (0+0w) 


(1) 


where  V 


B/Wt 


is  the  velocity  of  the  blast  relative  to  the  body.  In 


regions  1,  2,  and  3 of  Figure  6,  the  flow  velocity  components  in  the 
x'  and  y'  directions  relative  to  point  P are 

Ul/P  = Ul/W  " “P/W  U2/P  = U2/W  " “p/W  U3/P  = U3/W  “ “P/W 


'l/P 


1/W  ~ VP/W  V2/P  V2/W  VP/W 


= v„ 


- 


3/P 


V3/W  VP/W 


where 


“p/w 

= Vp/w  cos  ew 

vp/w 

= vp/w  sin  ew 

U1/W’  V1/W 

= Components  of 

U2/W’  V2/W 

= Components  of 

U3/W’  V3/W 

= Components  of 

_ . "I  f 

ub/vA 

relative  to 
the  body. 


and  6 = tan 


The  velocity  components  define  the  flow  direction  in  each  region. 
By  postulating  a pressure  jump  across  the  outer  segment  of  the  trans- 
mitted wave,  the  oblique  shock  relations  (e.g.,  see  Reference  11) 
permit  calculation  of  the  flow  direction  in  region  4.  The  same  shock 
relations  are  used  to  calculate  the  flow  direction  in  region  5,  on 
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the  condition  that  = P^.  This  iterative  procedure  is  continued 
until  the  flow  directions  in  regions  4 and  5 are  equal. 

The  outer  segment  of  the  transmitted  wave,  thus  determined,  is 
assumed  to  remain  straight  and  of  constant  strength  half-way  into 
the  shock  layer.  The  match  point  (R,  Figure  6)  location  and  velocity 
(Vr/W>  are  determined  by  the  point  of  intercept  between  the  outer 
segment  of  the  transmitted  wave  and  the  midpoint  of  the  shock  layer 
defined  by 

eM  - <ew  - ec)/2  (2) 

8i"  (VV 

and  VR/W  = Vp/W  sin  (E  -0  ) 

1 M 

Knowing  the  pre-interaction  conditions  at  the  center  of  the 
shock  layer,  the  pressure  behind  the  transmitted  wave  outer  segment 
is  defined.  The  compatible  inner  segment  is  found  by  assigning 
values  to  the  orientation  of  the  inner  segment  and  by  calculating 
the  resulting  pressure  jump  at  the  surface.  The  wave  strength  calcu- 
lated at  the  surface  is  presummed  to  apply  along  the  entire  inner 
segment  and  so  the  pressure  behind  the  inner  segment  at  the  center 
of  the  shock  layer  may  be  calculated.  The  solution  is  found  when 
the  pressure  calculated  behind  the  inner  segment  at  the  match  point 
agrees  with  the  value  calculated  behind  the  outer  segment.  In  this 
manner,  the  outer  portion  of  the  transmitted  wave  is  consistent  with 
the  blast-bow  interaction,  the  inner  segment  is  consistent  with  the 
wave-surface  interaction,  and  the  two  segments  are  compatible  with 
each  other  to  the  extent  that  the  pressure  behind  the  two  segments 
at  the  junction  point  is  matched. 

The  calculation  made  at  the  surface  begins  with  an  attempt  to 
find  a regularly  reflected  wave  which  satisfies  the  condition  that  the 
flow  must  be  parallel  to  the  surface.  If  none  is  found,  then  a Mach 
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stem,  which  is  infinitesimal  in  height,  is 
is,  the  stem  is  located  at  the  point  where 
transmitted  wave  impinges  on  the  surface, 
velocity  of  the  impingement  point  relative 


presumed  to  exist.  That 
the  inner  segment  of  the 
This  determines  the 

to  the  surface 

Q/W 


sin  (E  -0  ) 

V = V — — 

Q/W  R/W  sin  (Eo-0  ) 
L C 


(4) 


The  velocity  of  the  point,  Q,  relative  to  the  gas  on  the  cone  surface 
is  therefore  also  defined  which,  in  turn,  determines  the  strength  of 
the  wave  travelling  along  the  surface. 

The  Rankine-Hugoniot  equations  are  then  used  to  calculate  the 
pressure  jump  across  the  normal  shock.  The  key  assumptions  are  that 
the  transmitted  wave  is  not  so  curved  that  a two  piece  representation 
is  not  adequate  and  that  the  Mach  stem,  when  it  occurs,  is  small  in 
extent  relative  to  the  height  of  the  shock  layer.  In  principle,  the 
matching  of  pressure  only  at  the  point  of  bending  of  the  transmitted 
wave  is  incomplete  because  the  flow  directions  should  also  be  matched. 

It  was  not  possible  to  retain  the  simplicity  of  the  approach  and 
match  flow  directions  too.  Meanwhile,  the  comparisons  with  the  Ames 
peak  windward  pressures  at  a variety  of  conditions  show  agreement 
to  within  about  10  percent.  Furthermore,  a considerable  advantage  is 
gained  in  computation  time  because  the  new  model  requires  just  a few 
seconds  of  CPU  time  on  the  IBM  370/158  at  Martin  Marietta's  Orlando 
Division. 

With  two  exceptions,  the  new  model  treats  the  peak  windward 
pressure  from  nose-on  to  broadside  since  provision  was  made  for  the 
possibility  for  the  forward  moving  intercept  point. 

One  exception  occurs  when  the  blast  wave  is  exactly  parallel  to 
the  bow  wave.  Solving  the  interaction  between  the  blast  and  bow 
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waves  involves  defining  the  local  flow  velocities  relative  to  the  inter- 
section of  the  two  waves.  When  the  waves  are  parallel,  the  inter- 
section point  velocity  approaches  infinity,  so  a solution  is  not 
possible.  Loss  of  the  solution  at  the  single  value  of  encounter 
angle  is  not  seen  to  be  a serious  matter.  The  second  limitation  occurs 
at  the  transition  encounter  angle,  i.e.,  where  the  wave  pattern  changes 
from  Mach  stem  to  regular  reflections,  because  the  pressure  function 
appears  to  be  discontinuous.  Even  so,  the  error  in  pressure  is  on 
the  order  of  10  percent.  As  with  the  Ames  model,  the  new  model  has 
not  been  exercised  to  the  point  where  the  limits  of  applicability 
are  well  defined;  therefore,  a complete  judgment  of  this  point  cannot 
be  made  at  this  time. 

The  foregoing  analysis  was  computerized  and  used  to  calculate  the 
peak  windward  ray  pressures  at  the  set  of  conditions  which  were  tested 
during  the  Holloman  AFB  sled  test  program  described  in  Reference  2 (see 
Table  1) . For  a number  of  the  conditions,  predictions  are  also  calculated 
with  the  method  of  Reference  1 and  provided  by  Kutler  in  Reference  12. 


Table  1.  Sled  Test  Conditions,  Phases  II  and  III 


Test 

Number 

Cone 

Half-Angle(Deg) 

Mach 

Number 

Angle  of 
Attack  (Deg) 

Blast  Strength 
P2/P1 

Encounter  Angles  (Deg) 

B-2* 

11.2 

5 

0 

1.8 

22.2,  26.3  33,  78.4 

(Model  1) 

B-3 

2.2 

24.9,  28.8,  33.2,  78.1 

B-4 

2.2 

24.9,  31.8,  41.2,  78.2 

B-5 

2.6 

22.2,  28.9,  36.9,  86.5 

B-6 

24.7,  31.8,  45.0,  64.0 

B-7 

6 

24.8,  36.9,  44.9,  62.8 

(Model  2) 

L 

B-8 

11.2 

2.6 

24.8,  36.7,  45.0,  66.3 

(Mod 

el  1) 

B-9 

5 

5 

17.9,  21.9,  31.3,  78 

B-10 

■ 

■ 

5 

5 

** 

28,  32,  41 

*Test  B-l  was  aborted. 
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Figure  7 shows  the  peak  windward  pressures  plotted  against 
encounter  angle  (8)  Including  test  data,  Kutler's  results,  and  the 
curves  produced  with  the  model  discussed  here. 

A number  of  observations  are  possible  regarding  these  comparisons. 
Good  agreement  is  indicated  between  the  simple  model  and  the  rigorous 
solution  of  Kutler.  There  is  apparently  poor  agreement  between  the 
calculated  values  and  the  test  data.  However,  as  explained  in 
References  2 and  3,  the  test  data  are  expected  to  be  lower  than  the 
predictions,  particularly  in  the  Mach  stem  regime  because  the  pressure 
data  acquisition  system  used  in  the  sled  tests  did  not  have  a rapid 
enough  response  time  to  record  the  peak.  Reference  2 shows  that  when 
the  response  lag  is  considered,  the  Kutler  results  agree  very  well 
with  the  test  data.  Therefore,  the  Kutler  solution  is  treated  here 
as  the  standard  for  comparison. 

3.2  Peak  Leeward  Pressure 

Along  the  leeward  ray  of  the  cone,  the  transmitted  wave  pattern 
is  a Mach  stem  regardless  of  encounter  angle.  Furthermore,  the  Kutler 
results  show  that  as  the  encounter  angle  varies  from  nose-on  to  broad- 
side, the  peak  pressure  does  not  always  exceed  the  quasi-steady  value. 
Thus,  there  is  not  necessarily  a pressure  overshoot.  In  this  case, 
the  pressure  immediately  behind  the  transmitted  wave  is  less  than  the 
quasi-steady  pressure  and  it  rises  monotonically  to  the  quasi-steady 
value.  Therefore,  the  approach  taken  in  formulating  the  current 
leeward  pressure  model  is  to  combine  a method  for  calculating  the 
pressure  behind  the  transmitted  wave  with  a method  for  calculating 
the  quasi-steady  pressure.  Both  quantities  are  calculated  at  each 
encounter  angle  and  the  applicable  peak  pressure  is  the  higher  of  the 
two.  Typically  from  8=0  (nose-on  encounter)  to  some  6=B  , the 

pressure  behind  the  transmitted  wave  exceeds  the  quasi-steady 
pressure.  For  8>6qS»  the  peak  pressure  is  the  quasi-steady  value 


r 


T — 


f. 

T* 


and,  for  future  reference,  it  is  noted  that  the  pulse  duration  is 
therefore  zero.  As  discussed  earlier,  the  time  interval  in  question 
refers  to  the  time  it  takes  for  the  pressure  to  decay  from  the 
peak  to  the  quasi-steady  level.  Therefore,  there  is  no  shock-on- 
shock  pulse  over  a portion  of  the  leeside  of  the  cone  for  a portion 
of  the  encounter  angle  range,  and  that  portion  is  determined  by 
comparing  the  pressure  behind  the  transmitted  wave  to  the  leeward 
quasi-steady  pressure. 

The  pressure  behind  the  transmitted  wave  on  the  leeward  ray  is 
calculated  with  the  wave  pattern  approximated  as  shown  in  Figure  8. 
This  representation  is  simpler  than  that  used  on  the  windward  ray 
and  was  expected  to  suffice  because  there  is  no  transition  from 
Mach  to  regular  reflection  and  because  of  the  knowledge  gained  from 
the  windward  ray  analysis  at  8=0. 
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The  leeward  peak  pressures  calculated  with  the  simpler  wave 
approximation  are  used  to  determine  the  trend  with  encounter  angle  as 
B increases  from  zero.  Non-dimensionalizing  the  values  thus  obtained 
with  the  value  at  B=0  and  then  using  the  value  calculated  at  B=0 
with  the  windward  analysis  of  Section  3.1  produces  a good  estimate 
of  the  leeward  pressure  behind  the  Mach  stem  very  conveniently.  The 
leeward  analysis  proceeds  as  follows. 


The  pressure  behind  the  transmitted  wave  is  calculated  by  again 
using  the  Rankine-Hugoniot  relations  for  a normal  shock.  The  assump- 
tion that  the  transmitted  wave  is  straight  defines  the  position  and, 

therefore,  the  velocity  (V./tI)T  of  the  point  of  impingement  relative 

c^/w  Lee 

to  the  body  surface  as  given  by  the  following  expression. 


^VQ/vpLee 


V 


cos  (B~6) 

B/W  cos  (B-9  ) 
c 


(5) 


The  velocity  of  the  Mach  stem  relative  to  the  gas  on  the  cone 
surface  may  then  be  calculated  directly,  thereby  defining  the 
pressure  behind  the  stem. 


The  quasi-steady  pressure  (P  ) is  calculated  with  the  follow- 


ing  semi-empirical  expression  which  was  taken  from  Reference  13. 


QS 


P2  + J P2  V2  COS  Y 


p-r 


) cos  <f)  + 


<4*1 . 

m:  j 


(6) 


where  the  quantities  P2,  p2>  V2»  M2  are  tlie  Pressure*  density, 
velocity,  and  Mach  number  respectively  of  the  gas  behind  the  blast 
wave.  The  angle,  y,  is  the  angle  between  the  velocity  vector 
and  the  outward  normal  to  the  cone  surface. 


Figure  9 shows  the  peak  leeward  pressure  plotted  against 
encounter  angle  for  several  test  cases:  B-2,  B-5  and  6,  B-9  and  10. 
Comparisons  with  the  Kutler  results  are  included,  but  no  test  data 
are  presented  because  of  the  poor  quality  of  the  test  data  in  that 
region.  Agreement  between  the  two  methods  is  seen  to  be  quite  good. 
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3.3  Windward  Pulse  Duration 

The  formulation  used  to  model  the  duration  of  the  pulse  is  taken 
from  an  approach  suggested  by  Smyrl  (Reference  (14).  The  element  of 
gas  which  is  at  the  tip  of  the  cone  when  the  blast  arrives,  is 
flagged  and  tracked  as  the  wave  propagates  along  the  cone.  The  shock 
wave  is  travelling  faster  than  the  disturbed  (flagged)  gas  and  so  the 
distance  between  them  increases  with  increasing  time.  The  flagged 
element  of  gas  is  the  source  of  a secondary  disturbance  which  is  radi- 
ated at  the  local  speed  of  sound  (Figure  10).  At  any  point  in  time, 
the  extent  to  which  the  second  disturbance  has  radiated,  which  is 
defined  by  the  sonic  circle,  therefore  bounds  the  region  affected  by 
the  propagation  shock  wave. 


Figure  10.  A Model  of  the  Windward 
Pulse  Duration  - Regular  Reflection 

Smyrl  suggested  that  the  time  it  takes  for  the  pressure  to  reach 
the  quasi-steady  level  can  be  approximated  by  the  difference  in  time 
between  the  passage  of  the  wave  front  and  the  arrival  of  the  sonic 
circle  (point  S,  Figure  10).  This  reasoning  was  used  and  found  to 
produce  excessively  large  values  in  comparison  to  either  the  Kutler 
results  or  the  sled  test  data.  Therefore,  the  approach  was  modified 
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by  assuming  that  the  secondary  disturbance,  namely  the  one  emanating 
from  the  flagged  gas  element,  is  insignificant  and  that  the  quasi- 
steady pressure  is  reached  when  the  flagged  element  passes  the  point 
in  question.  This  approach  gives  much  better  agreement  with  the  other 
data  and  therefore  is  implemented  in  the  model  when  regular  reflection 
occurs. 

In  the  Mach  reflection  regime,  this  approach  was  found  to  give 
contradictory  results  in  that  the  pulse  duration  appeared  to  decrease 
with  increasing  £5.  The  Mach  stem  produces  a stronger  disturbance  and, 
thus,  the  secondary  disturbance  is  probably  not  insignificant  although 
the  duration  is  still  overpredicted  if  the  sonic  circle  reasoning  is 
used.  Therefore,  a new  model  was  formulated  as  follows. 

The  ratio  of  the  quasi-steady  pressure  to  the  pre-encounter 
surface  pressure  defines  an  effective  wave  front  which  follows  the 
actual  wave  front  (Figure  11).  The  speed  of  the  effective  wave  as- 
sociated with  the  quasi-steady  pressure  level  is  defined  by  this 
pressure  ratio  and  therefore  defines  an  estimate  of  the  arrival  of 
the  quasi-steady  level  of  pressure.  The  pulse  duration  is  approximated 
by  the  interval  between  the  arrival  of  the  actual  wave  front,  which 
is  known  from  the  peak  pressure  analysis  of  Section  3.1,  and  the 
arrival  of  the  effective  wave  front  associated  with  the  quasi-steady 
pressure. 

Figure  12  presents  comparisons  between  the  current  results, 
Kutler's  calculations,  and  the  sled  test  data.  The  figure  shows  the 
windward  pulse  duration  plotted  against  encounter  angles  for  the 
conditions  corresponding  to  sled  tests  B-2,  B-5  and  6,  B-9  and  10 
respectively.  It  should  be  noted  that  the  Kutler  and  sled  test  pulse 
durations  are  taken  from  pressure  time  traces  which  approach  the 
quasi-steady  value  asymptotically,  and  so  the  point  at  which  the  pulse 
ends  is  somewhat  subjective  and  therefore  is  not  always  precisely 
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defined.  However,  the  values  shown  are  representative  and  con- 
sistently arrived  at  and  do  lend  credibility  to  this  simple  model. 


re*K 


Figure  11.  A Model  of  the  Windward 
Pulse  Duration  - Mach  Reflection 

3.4  Leeward  Pulse  Duration 

A very  simple  model  was  chosen  for  the  leeward  pulse  duration 

because  of  the  information  available  from  the  analyses  previously 

performed.  It  was  noted  earlier  that  the  peak  leeward  pressure  at 

me  — encounter  angles  will  not  exceed  the  quasi-steady  pressure 

a therefore,  the  pulse  duration  is  zero.  The  encounter  angle 

at  which  the  crossover  occurs,  labelled  earlier  as  6qS,  is  typically 

no  larger  than  about  15  degrees.  That  is,  the  leeward  pulse  duration 

is  non- zero  only  between  6=0  (nose  on)  and  6=6n<3.  The  pulse  duration 

Wo 

at  6=0  is  already  defined  by  the  windward  analysis  and  so  the  leeward 

pulse  duration  is  simply  assumed  to  vary  linearly  with  6 from  the  value 

calculated  at  6=0  to  zero  at  6=6nc. 

QS 
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This  approach  is  illustrated  in  Figure  13  wherein  the  durations 
are  plotted  against  encounter  angle  for  the  conditions  corresponding 
to  sled  tests  B-2,  B-5  and  6,  B-9  and  10.  As  with  the  pressures 
on  the  leeslde,  there  are  no  sled  test  data  available  and  the  Kutler 

results  are  all  at  encounter  angles  beyond  3=8ric,  and  so  agreement  at 
8<8qS  cannot  be  checked.  However,  over  most  of  the  range  of  encounter 
angles,  where  the  duration  is  zero,  agreement  is  established. 

3.5  Circumferential  Distribution  - Peak  Pressures 

The  model  of  the  circumferential  distribution  of  peak  pressure 
uses  the  information  gained  from  the  windward  and  leeward  analyses 
regarding  the  position  and  velocity  of  the  wave  front.  The  assumption 
is  made  that  the  plane  defined  by  the  windward  and  leeward  wave  front 
locations  (see  Figure  14)  contains  the  wave  front  at  all  locations 
around  the  circumference.  The  line  of  intercept  between  the  plane, 
thus  defined,  and  the  cone  surface  is  taken  to  be  the  location  of  the 
wave  front.  The  wave  front  location  also  defines  the  wave  velocity 
relative  to  the  surface  (Vq)  as  given  by: 


V 


Q 


where 


K 


1 


(1  - K ) 

— (V  , ) 

cos  (j>  - K1  v Q/W'wind 


1 + ^VQ/W)lee/(VQ/W)  wind 
1 “ ^Q/W^lee^^Q/W^wind 


(7) 


(8) 


^VQ/W)wind’  ^Q/W^ee  = W3Ve  fr°nt  velocities  relative  to  the 

surface  along  the  windward  and  leeward 
rays  respectively. 

<f>  = circumferential  angle  measured  in  a plane  perpendicular 

to  the  cone  axis  with  <f>=0  at  the  windward  ray. 
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Figure  13.  Comparison  Between  the  New  Model  and 
Existing  Leeward  Pulse  Duration  Data 


The  wave  front  velocity  relative  to  the  gas  on  the  surface  in 
front  of  the  wave  may  then  be  calculated  which,  in  turn,  determines 
the  wave  Mach  number  and  pressure  jump  by  virtue  of  the  Raskine- 
Hugoniot  equations.  Since  the  angle  of  attack  is  restricted  to  small 
values,  the  effect  of  angle  of  attack  on  the  local  flow  direction  is 
neglected  in  computing  the  velocity  of  the  wave  relative  to  the  gas. 

It  is  recognized  that  the  actual  wave  front  will  actually  curve 
around  the  surface  so  it  will  not  be  in  a plane  as  suggested  here. 
However,  this  model  is  seen  to  produce  a plausible  distribution  of 
peak  pressure  as  evidenced  in  Figure  15  where  the  calculated  values 
are  plotted  against  the  circumferential  angle,  <j>.  The  values  obtained 
by  Kutler  and  those  measured  during  the  sled  tests  are  also  shown. 
Again,  the  comparisons  are  presented  for  the  B-2,  B-5  and  6,  B-9  and 
10  test  conditions,  and  reasonable  agreement  is  achieved  between 
Kutler' s results  and  the  current  model. 
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Figure  15.  Comparisons  Between  the  New  Model  and  Existing 
Circumferential  Pressure  Data 


3.6  Circumferential  Distribution  - Pulse  Duration 


The  approach  used  to  model  the  windward  and  leeward  pulse  durations 
is  not  easily  applied  to  the  circumferential  distribution  because  of  the 
difficulty  in  modeling  the  surface  streamlines  of  the  flowfield  behind 
the  blast  wave.  On  the  windward  and  leeward  rays  of  the  cone,  the 
streamlines  are  defined  because  they  coincide  with  the  geometric 
rays  of  the  cone.  At  other  circumferential  locations,  the  stream- 
lines are  curved  not  only  because  of  the  initial  angle  of  attack, 
but  also  because  of  the  angle  of  attack  induced  by  the  blast  wave. 

Thus,  attempting  to  follow  a flagged  element  of  gas  cannot  be 
achieved  without  solving  for  the  entire  flow  field  because  the  path  is 
not  known. 

An  approach  is  taken  which  is  analogous  to  that  used  to  model 
the  circumferential  distribution  of  peak  pressures  (see  Section  3.5). 

It  is  recalled  that  the  location  of  the  wave  front  around  the  cir- 
cumference is  approximated  by  the  surface  intercept  of  the  plane 
defined  by  the  windward  and  leeward  wave  front  locations.  In  a 
similar  manner,  the  windward  and  leeward  locations  of  the  quasi- 
steady region  define  a plane.  The  surface  intercept  of  that  plane 
around  the  circumference  defines  an  approximation  to  the  circumferen- 
tial location  of  quasi-steady  region.  The  pulse  duration  at  any  cir- 
cumferential location,  4>,  is  defined  as  the  time  interval  between 
the  passage  of  the  wave  front  and  the  passage  of  the  quasi-steady 
"front." 

Figure  16  illustrates  the  use  of  the  above  reasoning  in  sled 
tests  B-2,  B-5  and  6,  B-9  and  10.  Pulse  durations  are  plotted 
against  the  circumferential  angle,  <J>.  Unfortunately,  no  other 
results  are  available  for  comparison  from  either  the  Kutler  solution 
or  the  sled  tests. 
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Figure  16.  Typical  Distributions  of  Pulse  Duration 


3.7  Engulfment  Time 


An  additional  feature  of  the  phenomenon  which  is  of  interest 
to  the  structural  dynamist  is  defined  as  the  engulfment  time.  This 
quantity  pertains  to  the  passage  of  the  wave  front  over  a particular 
axial  station  of  the  body.  In  general,  the  wave  front  reaches  the 
subject  station  first  on  the  windward  ray  and  last  on  the  leeward 
ray.  The  differential  time  between  the  first  and  last  arrival  is 
called  the  engulfment  time  and  it  influences  the  manner  in  which  the 
energy  is  propagated  into  the  structure.  Obviously,  engulfment  time 
is  easily  calculated  from  the  results  produced  with  the  foregoing 
analysis  because  the  speed  of  the  wave  front  was  defined  at  the 
windward  and  leeward  rays.  Therefore,  the  engulfment  time  (ENGTIM) 
at  axial  station  X is  given  by 


(E) 


106X  f 


TIM 


cos  0r 


L(VQ/w\ee 


^Q/W^wind 


(9) 
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where 


%/W^lee 


cone  half-angle 

wave  front  velocity  relative  to  the  cone  along 
the  leeward  ray. 


^Q/W^wind 


wave  front  velocity  relative  to  the  cone  along 
the  windward  ray. 
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4.0  COMPUTER  PROGRAM 


The  analysis  described  in  the  foregoing  sections  was  programmed 
on  Martin  Marietta's  IBM  370/158  computer.  The  analysis  was  programmed 
using  FORTRAN  IV  through  a time  sharing  terminal  system.  The  following 
paragraphs  describe  the  structure,  input/output  features,  a complete 
listing  of  the  program,  and  a sample  problem. 

4.1  Program  Logic  and  Structure 

The  program  is  designed  to  analyze  the  encounter  between  a super- 
sonic sharp  cone  and  a planar  blast  wave.  The  cone  may  be  at  a small 
angle  of  attack  (a)  where  small  is  taken  to  mean  less  than  the  cone 
half-angle,  and  the  encounter  angle  (S)  may  be  anywhere  from  nose-on 
to  broadside.  The  only  encounter  in  this  range  which  cannot  be  treated 
is  the  case  where  the  blast  wave  is  parallel  to  the  bow  wave.  This 
is  not  a serious  matter  because  the  solution  is  well  behaved  at  smaller 
and  larger  angles  and  is  relatively  insensitive  to  encounter  angles  in 
that  regime. 

The  program  is  set  up  to  automatically  cover  the  entire  range  of 
encounter  angles  from  nose-on  to  broadside  with  a predetermined  incre- 
ment of  5 degrees.  This  arrangement  was  chosen  because  the  angle  at 
which  the  peak  pressure  is  highest  cannot  be  determined  a priori.  Since 
the  cost  of  running  the  program  is  small,  the  user  is  therefore  en- 
couraged to  allow  the  program  to  find  the  worst  condition  by  covering 
the  entire  range.  Since  there  is  a discontinuity  in  the  peak  pressure 
at  the  encounter  angle,  where  transition  from  regular  to  Mach  reflection 
occurs,  the  program  then  proceeds  to  perform  the  calculation  at  1-degree 
intervals  in  the  vicinity  of  the  discontinuity  to  reduce  the  likeli- 
hood of  missing  the  peak. 
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At  each  encounter  angle,  the  windward  peak  pressure  and  pulse 
duration  are  calculated  first  followed  by  the  analogous  calculations  on 
the  leeward  ray  of  the  cone.  The  circumferential  distributions  are  then 
evaluated  and  printed  out. 

The  analysis  is  performed  at  an  axial  station  (X)  selected  by  the 
user.  The  station  does  not  affect  the  pressure  calculation  because  the 
flow  field  is  presumed  to  be  independent  of  physical  scale  and  that  the 
wave  pattern  remains  self-similar  as  it  moves  down  the  body.  The 
pulse  duration,  on  the  other  hand,  does  vary  with  the  scale  of  the 
problem  because  it  is  determined  by  the  time  required  to  traverse  a 
physical  distance. 

Figure  17  is  a flow  chart  for  the  total  program.  The  initial 
decision  point  pertains  to  the  encounter  angle  range  of  interest 
followed  by  a choice  of  branches  for  zero  and  non-zero  angle  of  attack. 
Three  major  paths  are  identified  by  the  index  K.  The  first  (K=l) 
treats  the  case  of  6=0  for  the  leeward  ray  of  a cone  at  angle  of 
attack.  The  cone  at  angle  of  attack  is  treated  as  an  effective  cone 
at  zero  a with  cone  angles  equal  to  0C  + a and  6C  - a,  respectively, 
at  the  windward  and  leeward  rays.  At  6=0,  using  the  leeward  ray  flow 
conditions  and  the  windward  ray  analysis,  the  peak  pressure  and  pulse 
duration  are  calculated  and  stored  for  later  use  in  the  leeward 
analysis  at  6^0.  If  ot=0,  the  path  K=1  is  omitted. 

The  path  K=2  treats  the  complete  encounter  angle  range  from 
nose-on  to  broadside.  The  sequence  begins  with  the  analysis  of  the 
windward  ray  followed  by  the  leeward  ray  and  finally  the  circum- 
ferential distribution.  The  sequence  is  initialized  with  the  local 
flow  conditions  on  the  real  or  effective  cone  depending  on  whether  or 
not  the  angle  of  attack  is  zero. 

The  path  K=3  treats  four  additional  encounter  angles  in  the 
vicinity  of  the  value  of  6 at  which  transition  from  Mach  to  regular 
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Input /Output 
Input 

The  following  twenty  quantities  must  be  input  to  run  the  program 


Program  Symbol 


Definition 


N 

PI 

R1 

VI 

TC 

ALD 

P2P1 

TWOD 

PWP10 

RWR10 

V3WO 

TWLEED 

PWP1LE 

RWR1LE 

V3WLEE 


Number  of  blast  strengths 

2 

Freestream  pressure,  ///in 

3 

Freestream  density,  ///ft 

Freestream  velocity,  ft/s 

Cone  semi-vertex  angle,  deg 

Angle  of  attack,  deg 

Blast  pressure  ratio,  non-dimensional 

Bow  wave  angle,  ot=0,  deg 

Cone  surface-to-freestream  pressure  ratio 
a^O,  non-dimensional 

Cone  surface-to-freestream  density  ratio, 
a=0,  non-dimensional 

Flow  velocity  at  cone  surface  relative  to 
the  cone,  a=0,  ft/s 

Leeward  ray  bow  wave  angle,  a?/0,  measured 
from  the  effective  cone  centerline 

(Figure  4) , deg 

Cone  surface-to-freestream  pressure  ratio 
on  the  leeward  ray,  a^O,  non-dimensional 

Cone  surface-to-freestream  density  ratio 
on  the  leeward  ray  of  the  cone,  a^O; 
non-dimensional 

Flow  velocity  on  the  leeward  ray  of  cone 
ray  surface  relative  to  the  cone,  a^O, 
ft/s 
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TWWD 

Windward  ray  bow  wave  angle,  a^O,  measured 
from  the  effective  cone  centerline  (Figure 
15),  deg 

PWP1W 

Cone  surface-to-freestream  pressure  ratio 
on  the  windward  ray,  a^O,  non-dimensional 

RWR1W 

Cone  surface-to-freestream  density  ratio 
on  the  leeward  ray,  a^O,  non-dimensional 

VW 

Flow  velocity  on  the  windward  ray  of  cone 
surface  relative  to  the  cone,  ft/s 

X 

Axial  station  at  which  the  analysis  is 
performed,  ft. 

The  input  format  for  the  above  quantities  for  the  usual  batch 
processing  on  a computer  is  as  follows.  The  sequence  of  inputs,  as 
presented,  must  be  adhered  to  but  there  is  no  other  format  restriction 
other  than  the  total  number  of  quantities  on  each  card. 

Card  Input  Quantities 

1 N,  PI,  Rl,  VI,  TC,  ALD,  P2P1 

2 TWOD,  PWP10,  RWR10,  V3W0,  TWLEED , PWP1LE,  RWR1LE 

3 V3WLEE,  TWWD,  PWP1W,  RWR1W,  VW,  X 

Output 

The  output  information  (an  example  of  which  is  presented  in 
Section  4.4)  summarizes  some  of  the  input  data  and  contains  the  cir- 
cumferential distributions  of  peak  pressure  and  pulse  duration  at  each 
of  the  encounter  angles. 

The  input  quantities  printed  out  are  the  cone  semi-vertex  angle, 
freestream  velocity,  angle  of  attack,  axial  station,  and  blast 
pressure  ratio.  At  each  encounter  angle  (8)  that  the  analysis  is 
performed,  beginning  with  zero,  the  quantity  printed  first  is  the 
engulfment  time  (us)  followed  by  the  distribution  of  pre-encounter 
pressure,  peak  shock-on-shock  pressure,  quasi-steady  pressure,  and 
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pulse  duration  with  circumferential  body  angle  <p.  It  is  recalled 
that  <>=0  corresponds  to  the  windward  ray.  All  pressures  are  in  psi 
and  the  pulse  duration  is  given  in  ps. 

4.3  Program  Listing 

A listing  of  the  computerized  shock-on-shock  model  is  given  in 
the  following  pages.  For  the  benefit  of  the  user  and/or  programmer, 
a number  of  comments  are  included  in  the  listing.  The  comments 
spell  out  the  input  quantities  required  to  run  the  program,  to 
identify  major  subsections  of  the  analysis,  and  to  define  key 
variables. 

4.4  Sample  Problem 

A sample  problem  was  run  to  demonstrate  the  output  format  and  to 
provide  a check  case.  The  example  employs  the  conditions  correspond- 
ing to  sled  tests  B-9  and  B-10  as  listed  below. 


Quantity 

Value 

Quantity 

Value 

N 

1 

V3W0 

5426.0  ft/s 

PI 

12.6  psi 

TWLEED 

13.0  deg 

R1 

0.0627  lb/ft3 

PWP1LE 

1.69 

VI 

5500  ft/s 

RWR1LE 

1.396 

TC 

11.2  deg 

V3WLEE 

5584.0  ft/s 

ALD 

5.0  deg 

TWWD 

20.6  deg 

P2P1 

2.6 

PWP1W 

3.93 

TWOD 

16.6  deg 

RWR1W 

2.508 

PWP10 

2.596 

VW 

5375.0  ft/s 

RWR10 

1.95 

X 

1.25  ft 

To  reduce 

the  volume  of  paper,  the 

program  was 

modified  to  produce 

results  only  at 

four  encounter  angles: 

0,  20,  40, 

and  60  degrees. 

The  output  is  presented  on  the  facing  page. 
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SIMPLIFIED  S*CX-0*-SHDC*  MODEL 


CUNE  HALF- ANGLE-  I I .20 
AXIAL  S TATI  OH-  1.25 


CONE  VELOCITY-  5500,00 
BLAST  STRENGTH-  2,60 


A MOLE  OF  ATTACK- 


ENOULFMENT  TlMEl MICRO-SEC) - 


CIRCUMFERENTIAL 

ANGLE,PHI-UFO 

0.0 

20.00 

40.00 

60.00 

80.00 

100.00 

120.00 

140.00 

160.00 
I dO.OO 


PRE-FN COUNTER 
PRESSURE-PSI 

49.52 
4B.67 
4 6.22 

42.46 

37.86 

32.96 
2d.  35 

24.60 

22.15 
21  .29 


PEAK  PRESSURE 
PSI 

128.43 

126.29 
120.21 
I I I • 1 1 
100.28 
89.17 
79.  f | 

71  .17 
66.11 
64.37 


QUASI -STEADY 
PRESSURF-PSI 
117.98 
I 15.36 
108.12 
97.64 
85.85 
74.58 
65,20 
58.42 
54.41 
53.09 


FNGU LF -4 E NT  T IM F ( M I CRQ  - SE C ) - 


CIRCUMFERENTIAL 
ANGLE, PH I -DEO 
0.0 

20.00 

40.00 

60.00 

80.00 

100.00 

120.00 

I 40.00 

160.00 

180.00 


PRE-ENCOUNTFR 
PRESSURF-PSI 

49.52 
4d  .67 

46.22 

42.46 
37. d6 

32.96 
2d. 35 

24.60 

22.15 

21.29 


PEAK  PRESSURE 
PSI 

223.69 

215.26 
192.33 
160.83 
127.64 
98.20 
75.39 

59.87 
51  .08 

48.26 


QUASI -STEADY 
PRESSURE-PSI 
145.  J9 
140.42 
127.25 

108.74 

d8.76 

70.81 

57.05 

48.06 
43.27 
41  .81 


ENOULFMENT  T I *F.(  MICRO-SEC) » 


circumferential 

ANGLE,  Pill -!)F0 
0.0 
20.00 

40.00 

60.00 
dO.OO 

100.00 

120.00 

) 40.00 

160.00 
1 dO.OO 


PRE-FN COUNTER 
PRESSUPF-PSI 

49.52 

48.87 

46.22 

42.46 
37.86 

32.96 
28.35 

24.60 

22.15 
21  .29 


PEAK  PRESSURE 
PSI 

219.97 
2.0.29 

184.26 

149.29 
1 1 3.58 
83.11 
60.4b 
45.7| 
37.63 
35.09 


OUAS I -STEADY 
PRESSURE-PSI 

169.78 

162.79 
143.68 
117.33 

89.84 

66.45 

49.98 

40.5^ 

36.22 

35.09 


ENOULFMENT  TI4b< MlCRU-SEC)- 


PRc-ENaiUNTER 

PRESSURE-PSI 

49.52 

48.67 

46.22 

42.46 

37.86 

32.96 

28.35 

24.60 

22.15 

21 .29 


PEAK  PRESSURE 
PSI 

225.63 

215.23 
187.37 

150.23 
I 12.71 

81  .09 
57.94 
43.04 
34.98 

32.46 


OUAS I -STEADY 
PRESSURE-PSI 
186.64 
177.89 
154. JO 
121.77 
38.90 

62.13 

44.68 

35.99 

33.00 

32.46 
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5.0  CONCLUSIONS 


The  primary  conclusion  drawn  from  this  study  was  that  it  is 
possible  to  formulate  a reasonably  accurate  semi-empirical  model  of  the 
shock-on-shock  phenomenon.  As  with  any  such  model  the  range  of  appli- 
cability, until  it’s  demonstrated  otherwise,  is  assumed  to  be  limited  to 
the  range  of  parameters  for  which  test  data  are  available.  However,  the 
empiricism  was  kept  to  a level  which  should  allow  the  model  to  be  used 
effectively  over  a much  larger  range  of  conditions  (i.e.,  blast  strengths) 
than  that  covered  by  the  existing  data. 

The  model  is  cost  effective  in  that  a set  of  encounter  angles  can 
be  examined  with  no  more  than  a few  seconds  of  CPU  time  on  the  IBM 
370/158. 

The  model  deals  effectively  with  small  angles  of  attack,  but  must 
become  increasingly  suspect  as  the  angle  of  attack  exceeds  the  cone 
half-angle.  This  caution  is  noted  because  very  little  test  data  are 
available  for  a cone  at  angle  of  attack  in  a blast  environment  to 
compare  the  accuracy  of  the  prediction  method,  and  because  of  the 
occurrence  of  flow  separation  on  the  leeside  of  the  cone  when  the  angle 
of  attack  exceeds  the  cone  half-angle  (approximately) . The  occurrence 
of  separation  is  not  accounted  for  in  this  model  and,  therefore,  is  a 
possible  improvement  that  should  be  addressed  as  a future  task. 
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